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Abstract —Time domain characteristics of first and second 
order systems are well known. But the same simplicity and 
explicitness do not exist for low order fractional order 
systems (FOSs). Considering the step response, the 
templates are developed for designing the behavior of simple 
FOSs with a 2-term denominator polynomial (one is unity 
and the other involves fractional power). Although the 
explicit relations between design parameters and the 
performance parameters such as time constant, rise time, 
overshoot, settling time for fractional order control systems 
(FOCSs) do not exist and can’t be obtainable as in the 
ordinary integer order control systems, the obtained 
templates in this paper can be used for designing low order 
FOCSs. Hence, the drawback of non-existence of similar 
explicit formulas for FOCSs is eliminated by using these 
templates. 

Keywords — Control system , fractional order , rise time, 
step response, time constant. 

I. INTRODUCTION 

FOSs have founded many applications in the last two 
decades and a great deal of literature has appeared for 
analyzing and designing these systems [1-7]. Especially 
fractional order proportional integral derivative (FOPID) 
controllers have appeared very frequently in control system 
design [8-10]. 

Focusing on some very recent literature, for example [11] 
proposes an adaptive FOPI control method based on 
enhanced virtual reference feedback tuning to meet high 
precision and speed requirements for controlling flexible 
swing arm system in the light-emitting diode (LED) 
packaging industry. E. Cokmez et. all have obtained and 
visualized stability regions based on specified gain and 
phase margins for a FOPI controller to control integrating 
processes with time delay [12]. J. R. Nayak and B. Shaw 
have shown how to enhance the performance of the 
automatic generation control by adopting cascade 
proportional derivative (PD) - FOPID controller in a two- 


area mutually connected thermal power plant with 
generation rate constraint; group hunting search algorithm is 
adopted to explore the gain parameters of the controllers 
[13]. In [14], PI controller design is performed by using 
optimization for FOSs; first, controller parameters for a 
stable control are calculated by using the stability boundary 
locus method and then optimization is used to provide the 
best control. In [15], a new robust FOPID controller to 
stabilize a perturbed nonlinear chaotic system on one of its 
unstable fixed points is proposed based on the PID actions 
using the bifurcation diagram In [16], fractional-order 
discrete synchronization of a new fourth-order mends tor 
chaotic oscillator and the dynamic properties of the 
fractional-order discrete system are investigated; a new 
method for synchronizing is proposed and validated. 

In spite of the existence of a great deal of publications about 
FOSs some of which have just been mentioned above, most 
of the present analysis and design techniques deal with 
sophisticated and rather special applications [17-24]. 
Although the step response characteristics such as rise time, 
settling time, delay time, overshoot and some others are well 
known by explicit formulas for simple integer order systems 
[25], such formulas are not available for FOSs. And a 
compact publication yielding the relations between the design 
parameters and the step response characteristics of even simple 
FOSs are not yet present. The purpose of this paper is to 
fulfill this vacancy and to supply some design tools for 
simple order FOSs. 

The paper is organized as follows; in Section 2, basic 
definitions of time domain characteristics of first and second 
integer order systems are given. Section 3 introduces the 
FOSs that is studied and the investigation of its step 
responses depending on the fractional power. Section 4 gives 
and discusses the templates that can be used for the design of 
low order FOSs. Finally, Section 5 finishes with conclusions. 
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II. FIRST AND SECOND INTEGER ORDER 
SYSTEMS 

Let the first order system transfer function Hi be 

i 

tfiCs) = 


—■ (i) 

PlS + l 

It is assumed that p x > 0 for stability. Since the study is 
confined to step response characteristics, it is easily obtained 
by applying the unit step input u(t ) =0 for t < 0, u(t) = 


1 for t > 0, and compute the step response y(Y) as 

0, t < 0 


yW = [j. 


,-t/Pl 


t > 0 


( 2 ) 


The variation of the step response is shown in Fig. 1. It is an 
increasing exponential starting from 0 at t = 0, and rising to 
the steady-state value of 1 as lim t -> oo. The following time 
domain characteristics are defined for a response of the type 
shown in Fig. 1. 



Fig.l: Step response of a first order system and some important characteristics. 


Time constant t: It is the time required for the 
response to reach 1 - (l/e) = 0.632121 = 63 % of its 
final value. For this exponential, from Eq. (2) it is true that 

t = p 1 . (3a) 

Rise time T r : It is the time required for the response 
reach from 10 % to 90 % of its final value. 

T r — t 2 — ti = 2.302585 t - 0.105360r = rln9 

= 2.19722 5t = 2t. (3b) 

Settling time Ts: It is the time required for the 
response to stay around its final value with an errorless than 
2 %. 

T s = TlnSO = 3.912023t = 4t. (3c) 

Note that all the time characteristics depend on only the 
coefficient t = p ± in Eq. (1). So, p ± is chosen according to 
satisfy all the specifications on t ,T r ,T s . Note also that for 
p x = 0, the system is a unity gain system which yields 
y(t) = u(t)\ that is all the characteristic times z,T r ,T s are 
zero and no delay occurs at the response. 

In summary, the following properties of time domain 
characteristics are valid: i) The response increases 
exponentially to its steady state value without any 
oscillations; ii) Rise time and settling time are some 


multiples of time constant r = p 1 . 

It is well known that a time domain normalization with 
respect to T e corresponds to replacing t by t/r e which also 
corresponds to frequency domain normalization by writing s 
instead of sr e . In Eq. (1), if we replace s by s/p ± where 
T e = p x is the time constant, then we have 


ffiCs) = 


i 

5 + 1' 


(4) 


which has normalized time constant t = 1. 

As the reference transfer function for a second order system 

„> 2 

HXs) = 




5 2 + 2f (o n s + m 2 


( 5 ) 


is considered [25]; obviously dc gain is equal to 1; o) n is 
undamped natural frequency and > 0 is the damping ratio. 
For f > 1 which corresponds to overdamped case, Eq. (5) 
yields the step response 


1 / _J_ 

y(t) = 1- lr 2 e T 2 - r 1 e T i ), (6a) 

t 2 Ti V J 

which starts from 0 and rises to 1 monotonically as shown in 
Fig. 2 (zeta=2.0). In this expression, there are two time 
constants t 1 andr 2 so that 
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a> n O + y?^l)' 




(6b) 


"n(f-V? 2 -l)’ 

where si and si are the poles of transfer function. It can be 
shown by using (6b) that 


,_ <f 2^ 

+ t 2 = = -T n = —, (6c) 

7T U) n 

where T n = 2Z ) /oo n is the period of sustained oscillations of 
the undamped (f = 0) system. 



Fig.2: Step response of the second order system in Eqs. (6a), (6b) for values 2, 1, 0.2, 0 of damping ratio oo n = 1. 


For f = 1, the system in Eq. (5) is said to be critically 
damped; in this case there is only one time constant which is 
t = l/oo n = T n /2n. The step response is found as 


y(t) 


(Or, 


S 2 + 20)^5 + 00 2 s 


= 1 - + 1), (7) 

which increases monotonically from 0 to the steady-state 
value 1 as shown in Fig. 2 (zeta=1.0). 

For 0 < < 1, the system is said to be underdamped; for 

= 0.2 the step response is shown in Fig. 2 (zeta=0.2). It is 
seen that the response is stable and approaches to the 
reference value 1 in a damped oscillatory manner. It is a 

routine process to show that the response is given by 
e -Zo>nt 

y(t) = 1 - sin Ul~ % 2 (o n t 

V 1 - 

+ sm-yi^). (8) 


Finally, for <f = 0, the system is undamped and the step 
response is 

= = (9) 

which represents sustained oscillations as shown in Fig. 2 
(zeta=0) with undamped natural frequency oo n = 1 and 
undamped oscillation period T n = 2n. 


For the step responses corresponding to underdamped case 
two new time characteristics are defined. The oscillation 
period from Eq. (8) is 

2n n 

T 0 = = y (10a) 

V1 - f 2 ^n WW 2 

The settling time is approximately obtained from Eq. (8) as 
by equating the coefficient of sin function to 1 — 0.98 = 
0.02. The result is 

1 50 T e 50 

T s = -ln-= = -ln-^=. (10b) 

1 — % 2 2 f y]l — % 2 

Another time which is important is T max when the first peak 
occurs in the response. From Eq. (8) T max and y(T max ) = 
ymax are found t0 be 


T = 

1 max 


U 


n 


iZE 


ymax 1 


FT* 


(11a) 


(lib) 


Vw 2 

The overshoot y osh and percent overshoot (POSH) are 
defined by 


POSH 


yosh y-rr 


- 1 


yref 


■100 


100 e 


jw 1 


Vw 5 


(lie) 
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For the succeeding peaks, Eqs. (12a), (12b), (12c) are 
modified by replacing n by (2k — l)n where k represents 
the peak numbers; for the first peak k = 1. 

The reduction ratio (. RR ) is defined as the ratio of successive 
overshoots, and fromEq. (11c) 


2f7T 


RR = e ^ (lid) 

For the undamped system (( =0), Eq. (9) yields the 
following time domain characteristics: 


t = 1 + 


cos- 1 0 


OI„ 


(12a) 


cos UO.l) —cos 1 (0.9) 

T r = ---, (12b) 

0) n 

n 2n 

T max = —’ T s = 00 < T o = —• (12c, d, e) 

0)„ OJ n 

ymax 2, yosh 1' (12f, g) 

POSH = 100, RR = 1. (12h, i) 

For 0 < ( < oo, time domain characteristics 
T,T r ,T max ,T 0 ,T s and the overshoot y osh are plotted against 
f = r e /2.; The results are shown in Fig. 3. 



Fig.3: Variation of time domain characteristics against the damping ratio for % G [0,4-]. 


III. INVESTIGATION OF 2-TERM FRACTIONAL 
DENOMINATOR CHARACTERISTICS 

Consider the following fractional order transfer function 
with a constant numerator and 2-term fractional 
denominator: 

His) = - f . (13a) 

bs a + c 

Assuming dc gain (a/c) to be 1, letting b/c = p 1? and 


normalizing with r e 


. p u 

His) = 


1 we result with 
1 


s a + 1' 


(13b) 


We have the following observations on the step response of 
the fractional transfer function in Eq. (13b); see Fig. 4 for 
these observations: 
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Fig.4: Step response ofFOS in Eq. (14b) for different values of a. 


1) It is obvious that dc gain of this system is 1, which 
results from Eq. (13b) by inserting 5 = 0 and 
assuming a =£ 0. Therefore, all the step responses for 
stable cases will tend to 1 as limt^> oo; (a = 
0.5,1.2,1.7). 

2) For a = 0, the transfer function will be equal to a 
constant gain of 1/2, hence the step response is 0.5 
(<a = 0 ). 

3) For a = 1, the transfer function is equal to a first 
order integer type transfer function, hence, the step 
response is an increasing exponential with a time 
constant r = 1 (a = 1, see also Fig. 1). 

4) For a = 2, the fractional system in Eq. (13b) is 
equivalent to the second order integer type system in 
(5) with (o n = 1,^=0; hence, the step response is 
sustained oscillation (a = 2, see also Fig. 2). 

5) For a > 2, the system is not stable, and the step 
response increases exponentially (and os dilatory )-like 
manner (a = 2.018). 

6) For 1 < a < 2, the system is stable, and it has step 
responses (a = 1.2, a = 1.7). The first of these 
responses (a = 1.2) is a decaying curve after an 
overshoot; and the second (a = 1.7) is an oscillatory- 
like motion with exponentially-like decaying. 

7) For 0 < a < 1, step response is a stable 
exponentially-like increasing behavior (a = 0.5). 

We note that those responses for a = 2.018,1.7,0.5, 
resemble to those of a second order (for a = 2.018,1.7) and 
of a first order (a =0.5) integer order systems. But the 


explicit formulas as in Eqs. (3,10,11,12) between the system 
parameter a and the step response characteristics do not 
exist for the considered FOSs. So, in the following section, 
instead of using explicit formulas, some templates are 
obtained to be used for designing FOSs. 

IV. DEPENDENCE OF STEP RESPONSE 
CHARACTERISTICS ON a 

In this section the dependence curves (templates) of step 
response characteristics, namely, duration of first oscillation 
period (T 0 ), time constant t, rise time T r> and settling time T s 
vs a e [0.01,1.99], percent overshoot (POSH) vs /? = 2 - 
a E [0.01,1.99] are obtained by simulations. Simulations 
are carried for 30 5in steps of A a = 0.01 by subprograms of 
FOMCON toolbox [7] integrated with MATLAB R2017 [7]. 

Fig. 5 shows the variation of the duration of the first 
oscillation against a. Numerical data show that the first peak 
occurs for a = 1.01 and it is equal to 1. 0014. Then, until 
a = 1.34 second peak does not appear; more elaborate 
numerical analysis show that, the second maximum starts 
exhibiting for the first time for a = 1.3396 for which the 
first andsecond overshoots are 0.63695221,0.011386081, 
respectively; but for a =1.34 following the first peak of 
value 1.1640, the second peak of value 1.0114 occurs. This 
means period of the first oscillation is defined for a > 
1.3396. Since there are no peaks (maximums) until a = 
1.01, the graph is started from a = 0.8, though numerical 
data is obtained for all a £ [0.01,1.99]. 
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Fig. 6 shows plots of time constant t(Tclu ), rise time T r , 
and settling time T s . Time constant plot Tan starts from a = 
0.14, because, for smaller values of a,r is larger than 30 5 
so that the response can’t reach the critical value 0.632121 
until 30 s; similar arguments are true for the rise time T r 
which starts at a = 0.51, and for settling time T s which 
starts at a = 0.78. Tan is obviously decreasing with 
increasing a, and it becomes 1.1915 for a = 1.99. T r also 
decreases with increasing aand it changes from 28.0350 at 
a = 1.51 to 1.0220 at a = 1.99. For a = 1.70, t is almost 
equal to T r (r = 1.0965, T r = 1.0970); fora G [0.01,1.70], 
r < T r ; and for a E (1.70,1.99], t > T r . Settling time plot 


starts from T s = 27.33 for a = 0.78 and ends at T s = 29.82 
for a = 1.99. T s decreases until a = 1.07 and reaches to its 
minimum value T s = 2.8345 at a = 1.07, then it jumps up 
to T s = 5.1525 at a = 1.08. The plot terminates at a = 
1.99 with T s = 29.8215 . The irregular shape of increase of 
T s for a E [l.07,1.99] is due to the dependence of T s on 
discrete change of oscillations remaining in the 
limit[ 1, + 0.02]. Contrary to monotonic decrease of r and 
T s , it is true that (disregarding the irregular changes 
mentioned) T s decreases monotonically for a E [0.01,1.07] 
and it increases for a G [1.07,1.99], 



alpha 


Fig.6: Time constant t, rise time T r , and settling time T s vs a. 


Fig. 7 shows the comparison of the variation of overshoots 
with a G [1.01,1.99] (/? = 2 — a,/? G [0.01,0.99]) and 


with the damping ratio ^ of a second order system. This plot 
is useful for finding the fractional order a and damping ratio 
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( for a given overshoot. For example, to achieve an are appropriate for FOS and for a second order system, 

overshoot of 60 %; 2 — a = 0.24 -> a = 1.76 and = 0.16 respectively. 



Fig.7: Overshoots versus j3 = 2 — a (for fractional) and (for 2nd order) systems. 


Fig. 8 better illustrates the relations between the overshoot for a second order system, 

and directly a (not /? = 2 — a) for fractional system, and ^ 



Fig. 8: Variation of a for fractional order, and df for second order systems withthe overshoot. 


Fig. 9 shows the values (3 = 2 — a and against the rise 
time. It is obvious that for a rise time of 4.371, ft and ( have 


the same values of 1.2. This means for = 1.2 and fora 
2 — (3 = 2 — 1.2 = 0.8 the rise times are equal to 4.371. 
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Fig. 9: /? = 2 — a versus rise time and versus rise time. 


V. CONCLUSIONS 

Time domain characteristics of the FOS with a 2-term 
denominator polynomial involving a single fractional power 
is investigated in this presentation. Dependence of important 
step response characteristics, namely rise time, settling time, 
delay time, overshoot, and oscillation period on the 
fractional order a are derived, and the results are presented 
in graphical forms that can be used as templates for design 
purposes. The study is conducted comparatively by 
considering integer order systems of 1 st and 2 nd order types. 
It is shown that the same simplicity and exp licitness present 
for second order systems do not exist between the transfer 
function parameters and the step response characteristics for low 
order fractional systems. The results bring light for 
designing simple FOCSs, thus a vacancy has been fulfilled 
by this work. 
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